A0-A  C83S7/ 


A  J^UiXiN  xOxxJU 

t  library 


TECHNICAL  REPORT  ARLCB-TR-80002 


STRESS  SINGULARITY  AT  THE  VERTEX  OF  A  FLAT  WEDGE-SHAPED  CRACK 

.  BY  A  VARIATIONAL  METHOD 


M.  A.  Hussain 
B.  Noble 
S.  L.  Pu 


January  1980 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 

LARGE  CALIBER  WEAPON  SYSTEMS  LABORATORY 
BENI:T  WEAPONS  LABORATORY 
WATERVLIET,  N.  Y.  12189 


AMCMS  No.  611102H600211 
DA  Project  No.  1L1161102H60 
PRON  No.  1A924324GGGG 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


OTIC  QUALITY  INSPECTED  1 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other  author¬ 
ized  documents. 

The  use  of  trade  name(s)  and/or  manufacturer(s)  does  not  consti¬ 
tute  an  official  indorsement  or  approval. 


DISPOSITION 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it 


to  the  originator. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whin  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

t.  REPORT  NUMBER 

ARLCB-TR- 80002 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtltte) 

S.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

STRESS  SINGULARITY  AT  THE  VERTEX  OF  A  FLAT  WEDGE- 

SHAPED  CRACH  BY  VARIATIONAL  METHOD 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORO) 

M.  A.  Hussain 

B.  Noble 

S.  L.  Pu 

8.  CONTRACT  OR  GRANT  NUMBERffl) 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 

AREA  a  WORK  UNIT  NUMBERS 

Benet  Weapons  Laboratory 

AMCMS  No.  611102H600211 

Watervliet  Arsenal,  Watervliet,  NY  12189 

DA  Project  No.  1L1161102H60 

DRDAR-LCB-TL 

PR0N  No.  1A924324GGGG 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  REPORT  DATE 

US  Army  Armament  Research  and  Development  Command 

January, 1980 

Large  Caliber  Weapon  Systems  Laboratory 

13.  NUMBER  OF  PAGES 

Dover,  New  Jersey  07801 

30 

14.  MONITORING  AGENCY  NAME  ft  ADDRESS^//  different  from  Controlling  Office) 

IS.  SECURITY  CLASS.  (o{  this  report) 

UNCLASSIFIED 

1 DECLASSI  FI  CATION/ DOWN  GRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thla  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  abatract  entered  In  Block  20,  if  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 

Presented  at  1979  ASME-CSME  Applied  Mech/Bioengineering  Conf/Fluid  Eng  Conf, 
Niagara  Hilton  Hotel,  Niagara  Falls,  NY,  18-20  June  1979. 

Presented  at  2nd  MACSYMA  Users'  Conf,  Washington,  DC,  20-22  June  79. 

To  be  published  in  an  open  literature  journal. 


75!  KEY  WORDS  (Continue  on  revere*  etde  if  neceeeery  and  Identity  by  block  number) 

Variational  Method  Spherical  Harmonics 

Flat-Wedge-Shaped  Crack  Mixed  Boundary  Value  Problems 

Stress  Intensity  Factors  *  v  l 

Fracture  Mechanics  ~ 

Papkovitch  Stress  Functions 


20,  ABSTRACT  (VaaOmua  an  rereraa  aMt  Ft  nwMHqr  mod  Identify  by  block  number) 

Three  dimensional  elasticity  problems  are  generally  complex.  In  this  paper 
we  present  the  analysis  for  the  stress  singularity  at  the  apex  of  a  three  di¬ 
mensional,  flat,  wedge-shaped  crack  under  general  loadings.  The  problem  is 
reduced  to  a  set  of  coupled  dual  integral  equations.  Because  of  the  complex¬ 
ity  they  are  not  amenable  to  a  closed  form  solution.  A  variational  method  is 
developed  to  handle  such  problems.  The  physical  interpretation  of  the  results 
is  also  presented. 


DO  .SE.  H73  EDITION  OF  I  NOV  68  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Entered) 


SECURITY  CLASSIFICATION  OP  THIS  PAGEfWTi  an  Data  Entered) 


SECURITY  CLASSIFICATION  OF  THIS  PAGEf^en  Data  Entered) 


TABLE  OF  CONTENTS 

Page 

INTRODUCTION  1 

VARIATIONAL  THEOREM  2 

PAPKOVITCH  STRESS  FUNCTIONS  4 

THREE  MODES  OF  CRACKS  AND  COUPLED  INTEGRAL  EQUATIONS  6 

MODE  I  7 

MODE  II  8 

MODE  III  8 

APPLICATION  OF  VARIATIONAL  PRINCIPLE  11 

CONCLUSIONS  14 

REFERENCES  15 

APPENDIX  1  16 

APPENDIX  2  19 

TABLES 

I.  BOUNDARY  CONDITIONS  ON  0  =  ir/2  7 

II.  VALUES  OF  y  FOR  V  =  0.25  13 

ILLUSTRATIONS 

1.  A  flat  wedge-shaped  crack  under  three  different  modes  and  23 

its  two-dimensional  counterparts. 

2.  A  spherical  coordinate  surrounding  the  apex  of  a  thin  24 

angular  sector  crack. 

3.  The  eigenvalue  y  as  a  function  of  a/ir  for  v  =  0.25. 

(a)  Variational  results  indicated  by  x  and  •  for  mode  II.  25 

(b)  Variational  results  indicated  by  x  and  •  for  mode  III.  26 


i 


INTRODUCTION 


The  theory  of  fracture  mechanics  has  been  a  very  successful  tool  in 
engineering  application  in  recent  years.  This  is  mainly  due  to  the  use  of  a 
single  characteristic  parameter  namely  the  stress  intensity  factor,  that  is 
the  coefficient  of  the  stress  singularity  at  the  tip  of  a  crack  in  the  linear 
theory  of  elasticity,.  In  most  of  the  two  dimensional  cracks,  in  homogeneous 
media,  the  singularity  is  of  the  order  one  half.  For  the  three  dimensional 
cracks,  however,  the  singularity  depends  upon  the  geometric  configurations. 

In  this  report  we  study  the  singularities  at  the  apex  of  a  thin  wedge- 
shaped  crack  shown  in  Figure  1  under  three  loading  conditions.  Using  the  near 
field  approach,  the  problem  is  reduced  to  an  eigenvalue  problem  for  coupled 
dual  integral  equations.  The  results  indicate  that  cracks  tend  to  straighten 
out  at  the  apex. 

In  the  following  sections  we  first  prove  the  variational  theorem  by 
which  the  eigenvalue  of  coupled  integral  equations  is  to  be  obtained.  This 
eigenvalue  problem  is  not  the  linear  one  commonly  encountered  in  mathematical 
physics.  Next  we  present  Papkovitch  stress  function  approach  to  three  dimen¬ 
sional  theory  of  elasticity.  Then  coupled  dual  series  relations  are  obtained 
by  using  mixed  boundary  conditions.  These  series  are  transformed  to  coupled 
singular  integral  equations.  Finally  the  variational  method  is  applied  to  the 
coupled  integral  equations  to  obtain  eigenvalues. 
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The  variational  method  completely  avoids  the  solution  of  complex  singular 
integral  equations .  In  this  study  Macsyma  was  found  to  be  an  indispensable 
tool  at  all  levels  of  analysis. 


VARIATIONAL  THEOREM 

Consider  the  following  pair  of  homogeneous  coupled  integral  equations 
with  Fredholm  kernels. 

a  a 

/  K  C^JiOfOJO#  +  /  K12(<j>,ijsii)gO!W  =  0  (l) 

o  11  o 

a  a 

/  K21(<t,iP;U)f(ilj)dip  f  K22(4>,iJ;;v0g(iJ;)chJj  =  0  (2) 

0  0 

where  eigenvectors  f  and  g  and  the  eigenvalue  y  are  unknown  and  (4>,i|>;y) 

etc,  involve  y  in  a  linear  or  nonlinear  fashion*,  Construct  the  following 

characteristic  equation  for  the  determination  of  y*  with  appropriate  trial 

functions  f*(t|0  and  g*W* 

a  a  a  a 

(/  f*(x)/  K,,  (x,y) f* (y)dydx)  (/  g *(40/  K22M>,<i0g* 

0  0  0  0 

a  a  a  a 

-(/  f*(x)/  K  (x,y)g*(y)dydxK/  g*(<f0/  (<p,ip)  f*  (ip)dipd<p)  =  0  (3) 

0  0  0  0 

If  f*  and  g*  vary  around  the  exact  solutions  f  and  g  as 

«  f(40  +  65  M0  ,  =  gMO  +  6nM0  (4) 

then  (y*-y)  is  stationary  around  6  as  6  approaches  zero.  That  is 

y*  *  y  +  Offi2) 

provided 


2 


(5) 


KuOM)  =  kuGM) 

K22W,*)  =  K22(i|>,(J>) 

k12C<J>,^)  =  k21(<M) 

Proof:  Symbolically  we  write  equation  (3)  as 

(f*Knf*)(g*K22g*)  -  (f*K12g*)(g*K21f*)  =  0  (6) 

where 

f*  =  f  +  ,  g*  =  g  +  .  (7) 

Expanding  the  kernels  around  y,  we  have 

Kn(u)  =  Kn(y*)  +  AK«u(y*)  +  0(A2)  ,  etc.  (8) 

where 


A  =  y  -  y* 


(9) 


Substituting  from  (7)  and  (8)  into  (6)  and  using  (5)  we  obtain 


(CfKufMgK22g)  -  CfK12g)CfK12g)] 

♦  2«[(CKnf)(gK22g)  -  (5K12gHfK12g)  *  CnK22g)(fKnf)  -  (fK^rfl  (fK12g) ) 

+  L^A  +  L262  +  higher  order  terms  =  0  (10) 


where 


L1  =  (£KnO(gK'22g)  +  (£K»nf)(gK22g)  -  2(fK12g)(fK'12g)  (10a) 

L2  ■  (5KnC)  (gK22g)  -  (5K12g)2  *  4(fK11CMgK22rt 

*  (fKuf)(nK22n)  -  (ft12n)2  -  4(«12n) ta12g)  •  (i«b) 

Using  equations  (1)  and  (2)  it  is  seen  that  the  first  two  terms  in  (10)  vanish. 


Hence 

A  =  62(L2/L1)  +  0 ( 6 3 ) 


(ID 
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The  above  proves  that  £  is  of  the  second  order  in  60  This  completes  the 
proof  of  the  theorem.  Because  the  eigenvalues  involved  in  kernels  can  be  in 
a  nonlinear  fashion  we  cannot  prove  its  bounds  as  can  be  done  for  the  linear 
case  in  Reference  1. 


PAP KOV ITCH  STRESS  FUNCTIONS 

In  the  absence  of  body  force  the  equation  of  equilibrium  of  a  homogeneous 
isotropic  elastic  solid  is  given  by 

1  _  - 

(12) 


V2u  +  — —  VV*u  =  0 


l-2v 

Using  Helmholtz  decomposition  theorem  the  solution  of  (12)  can  be  written 
2,3 


as 


2Gu  =  VBq  +  V (r° B)  -  4(l-v)B 


(13) 


where  G  and  v  are  shear  modulus  and  Poisson's  ratio,  B  and  B  are  known  as 

o 

Papkovitch  stress  functions  satisfying 

V* 2Bq  =  0  ,  V2B  =0  ,  B  =  iip  +  jo)  +  kA  .  (14) 

For  computational  purposes  it  is  convenient  to  write  the  solution  as  a  super¬ 
position  of  the  following  basic  solutions. 

1st  Basic  solution:  2Gu  =  VBq 

2nd  Basic  solution:  2Gu  =  V(xip)  -  4(l-v)ijji 

3rd  Basic  solution:  2Gu  =  V(yw)  -  4(l-v)wj 

4th  Basic  solution:  2Gu  =  V(zA)  -  4(l-v)Ak 


(15) 


■*-Barlett,  C.C.  and  Noble,  B.,  "A  Variational  Method  for  the  Solution  of  Eigen¬ 
value  Problems  Involving  Mixed  Boundary  Conditions,"  Applied  Science  Research, 
Section  B,  Vol.  9,  1962. 

2Green,  A.E.  and  Zerna,  W. ,  THEORETICAL  ELASTICITY,  2nd  Edition,  Oxford,  1968. 

^ Lure ' ,  A.T.,  THREE  DIMENSIONAL  PROBLEMS  OF  THE  THEORY  OF  ELASTICITY,  Inter- 
Science  Publishers,  1964. 
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These  basic  solutions  were  transformed  to  the  following  spherical  coordinate 
system 

x  =  r  sin  0  cos  cj)  ,  y  =  r  sin  0  sin  <j>  ,  z  =  r  cos  0  (16) 

with  the  origin  at  the  apex  of  the  crack  as  shown  in  Figure  2.  This  enables 
us  to  obtain  near  field  solutions  and  to  study  the  stress  singularity  at  the 
apex.  The  complete  results  of  components  of  displacement  and  stress  for  basic 
solutions  are  given  in  Appendix  1.  Since  we  are  interested  in  the  power 
singularity  at  the  apex,  so  we  choose  Papkovitch  potentials  in  the  form 

H(r,e,40  =  r^Hj  (0,<t>)  (17] 

where  y  is  the  eigenvalue  to  be  determined.  As  can  be  seen  from  Appendix  1, 
stress  components  will  be  of  the  form 

a  =  0(ry_1]  (18) 

which  will  be  singular  when  y  <  1.  For  the  displacements  to  be  finite  we 
seek  positive  eigenvalues  between  0  and  1. 

The  near  field  geometry  surrounding  the  apex  permits  us  to  write  (17) 
as  the  separation  of  variables  solution 

r  u  m  cos  ^ 

H(r,e,4i)  =  l  rMP  (cos  6)  or  (19) 

m  s  in  m<j) 

where  P™(x)  is  the  associated  Legendre  function  of  the  first  kind  with  degree 
y  and  order  m.  When  (19)  is  substituted  into  the  four  basic  solutions  given 
in  Appendix  1  we  see  that  the  forms  of  some  solutions  thus  obtained  are  not 
convenient  to  work  with.  The  final  solutions  used  in  this  analysis  are 
designated  as  solutions  A,  B,  C,  and  D  which  are  given  in  Appendix  2.  Solu¬ 
tion  A  was  obtained  by  replacing  y  by  y  +  1  after  the  substitution  from  (19) 
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into  the  first  basic  solution.  Solution  B  is  obtained  by  adding  the  second 
and  the  third  basic  solutions  with  proper  trigonometric  functions  in  (19), 
replacing  m  by  in  +  1  and  then  using  Legendre  recursion  formulae.  A  similar 
process  was  used  to  obtain  solution  C.  Solution  D  simply  comes  from  the 
fourth  basic  solution  with  the  use  of  (19) 0 

Solutions  A,  B,  C,  and  D  must  not  be  linearly  independent.  This  is  due 
to  the  fact  that  the  condition  that  the  vector  point  function  is  solenoidal 
has  not  been  used  explicitly  in  Papkovitch  stress  functions  approach.  We 
found  these  solutions  are  indeed  not  independent.  A  relation  among  them  can 
be  written  symbolically  in  the  form 

(y+m+1)  [C]  =  [B]  +  2  (y+m+1)  [D]  -  2(y-3+4v)[A]  .  (20) 

Hereafter  solution  C  is  replaced  by  solutions  A,  B,  D  using  (20). 

THREE  MODES  OF  CRACKS  AND  COUPLED  INTEGRAL  EQUATIONS 

For  a  crack  shown  in  Figure  2,  the  leading  edges  of  the  crack  are  <J>  =  ±  a 
and  the  crack  is  in  the  x-y  plane  (0  =  tt/2)  0  Let  D~  and  D+  be  the  cracked  and 
uncracked  region  of  the  plane  9  =  7t/20  Within  the  cracked  region,  the  dis¬ 
placement  is  discontinuous.  If  the  discontinuity  is  in  the  z-direction 
(u*  “  u”  =  finite),  the  crack  is  under  mode  I;  if  the  discontinuity  is  in  the 

u  u . 

+  — 

x-direction  (u  -  u  =  finite),  the  crack  is  defined  to  be  under  mode  II;  and 
x  x 

if  u+  -  u  =  finite,  the  crack  is  defined  to  be  under  mode  III.  Boundary 

y  y 

conditions  for  various  modes  are  tabulated  below. 
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TABLE  I.  BOUNDARY  CONDITIONS  ON  9  =  tt/2 


Non  Mixed  Conditions  Mixed  Conditions 

fin  D  +  D  ) in  D~ in  D+ 


Mode  I 

TSr  '  TS4,  =  0 

o 

II 

CD 

c 

CD 

II 

o 

(21) 

Mode  II 

Q 

CD 

tl 

o 

tq  =  xQ .  =  0 

0r  64> 

u  -  u ,  =  0 

r  (J) 

(22) 

Mode  III 

Q 

CD 

II 

O 

T  =  T„  =  0 

0r  0(J> 

=  u .  =  0 

r  <p 

(23) 

For  mode  I,  Ug  is  even  in  <p.  This  leads  to  the  use  of  the  trigonometric 
functions  at  the  top  of  (19).  The  boundary  conditions  of  (22)  and  (23)  are 
identical,  but  the  symmetric  properties  are  different  for  mode  II  and  mode  III. 
In  the  former  case,  ur  is  even  and  u<j>  is  odd  in  <J>  while  in  the  latter,  the 
reverse  is  true.  Hence  the  proper  set  of  quantities  should  be  selected  in 
(19),  for  each  case,, 

MODE  I.  Using  equation  (19)  and  the  non-mixed  conditions  of  (21),  we 

have 

Bm  =  0  »  \  =  0n^l)'* 1(l-2v)Dm  .  (24) 

The  mixed  boundary  conditions  of  (21)  and  using  (24)  and  (19),  yield 

l  bm  cos  m<J>  =  0  0  <_  cj>  <  a 

r  (25) 

)  q„b  cos  md)  =  0  a  <  d>  <  tt 
u  m  T  r  — 

where  £  denotes  the  summation  with  respect  to  m  for  m  =  0,1,2,...°°,  and 

bm  =  (-m+y+l)DmP™+1  ,  qm  =  =  P°(0)  (26) 
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MODE  II.  For  the  homogeneous  condition  of  Og  in  (22),  using  (19),  the 
coefficients  A,  B,  and  D  are  related  in  the  form 

2(1_V)Dm  =  Cm+y+l)\  +  [y2  +  2v(y+l)  -  2(l-v)m  -  l]Bm 
This  relation  and  the  mixed  conditions  of  (22)  yield  the  following  coupled 
dual  series: 


I  Em  cos  m<t>  =  0 

T  (R  E  +S  F  )cos  mc|)  =  0 
L  m  m  m  m' 

r» 

l  F^  sin  m<()  =  0 

y  (U  F  +T  E  )sin  4=0 
u  m  m  m  nr 


0  <  (Ji  <  a 

a  <  <J>  <_  tt 

0  <_  $  <  a 

a  <  d>  <  it 


where  £  denotes  the  summation  with  respect  to  m  for  m  =  1,2,...“,  and 

E  =  - (m+y+1) {yA  +  [4m(l-v)2  +  y (y+4v-m-3) ] E  }P“ 
m  m  ^  J1  m  y 

F^  =  (m+y+1) {mA^  +  [4y(l-v)2  +  m(y+4v-m-3) ] B^JP™ 


(27) 


(28) 


(29) 

(30) 


R  =  [m2  -  y (y+1) (1-v) ] V  ,  S  =  m(l-v-vy)V 
m  m  m  m 

Um  =  [ (l-v)m2  -  y(y+l)]Vm  ,  Tm  =  mfl+vy)^ 

in  (31)  Vm  stands  £or  p™+i/py/ [ (rn+y+1) (m2-y2)] . 

MODE  III.  Similar  to  the  preceding  case,  we  have 

sin  m(j)  =  0  0  <  cf)  <  ct 

y  (R  E  +S  F  )sin  4  =  0  a  <  <b  <  tt 

m  m  m  m  — 

\  Fm  cos  mcf)  =  0  0  <  cj>  <  cx 

T  (U  F  +T  E  )cos  mcb  =  0  a  <  cb  <  7T 

L  v  m  m  m  nr  Y  Y  — 


(31) 


(32) 


(33) 
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It  was  found  that  the  dual  series  (25)  for  mode  I  is  identical  to  that 
of  a  potential  problem  studied  in  References  4,  5,  and  6  and  will  not  be 
discussed  here0 


For  convenience  we  make  the  following  change  of  variables 


(f)  =  7T-U3 


,  01)mE  =  E  » 

K  m  m 


(-l)mF  =  F  1 
m  m 


R  =  R  1 (1-v-vp) 
m  m  v  J 

u  =  u  1  (.l+vtOC-i) 

m  m  v  K  J 


Sm  =  S' (1-v-vy) 
m  m  ^  H J 

T  =  T  1 (l+vu) (-1) 
m  m  ^  K  J 


The  dual  series  of  (27) ,  (28)  now  become 

y  *  cos  moa  =  0 
L  m 

7  F  '  sin  mo)  =  0 
L  m 


a1  <  03  <  TT 


y  (R  1 E  *+S  1 F  !)cos  m03  =  0 
u  m  m  m  m  J 

y*  (T  1 E  1 +U  1 F  f)sin  mo)  =  0 
L  mm  mm 


0  <  a)  <  a1 


(34) 

(35) 

(36) 

(37) 

(38) 


Let  the  right  hand  sides  of  (35)  and  (36)  for  the  interval  0  <  03  <  a'  be 
denoted  by  unknown  functions  f(o3)  and  g(o3),  respectively.  The  Fourier 
inversion  gives 


a1 


ot 1 

Eq,  =  I  /  fOJO#  , 


Era’  =  7  /  f  010  cos 


0 


0 

a' 


(39) 


Fm'  =  £  J  g  010  sin 


^Noble,  B.,  "The  Potential  and  Charge  Distribution  Near  the  Tip  of  a  Flat 
Angular  Sector,"  EM-135,  New  York  University,  NY,  1959. 

^Brown,  S.N.,  and  Stewartson,  K.,  "Flow  Near  the  Apex  of  a  Plane  Delta  Wing," 
Journal  of  Institute  of  Mathematics  and  Its  Applications,  Vol.  5,  p.  206, 
1969. 

^Morrison,  J.A.,  and  Lewis,  J.A.,  "Charge  Singularity  at  the  Corner  of  a  Flat 
Plate,"  SIAM,  Journal  of  Applied  Mathematics,  Vol.  31,  p.  233,  1976. 
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Substituting  from  (39)  into  (37),  (38)  and  interchanging  the  order  of  summa¬ 
tion  and  integration  we  have  the  following  coupled  integral  equations  for 
the  determination  of  f  and  g. 


a'  a' 

/  K11(u,^;y)f (ip)d\p  +  /  Kl2  (aj,ip;y)g(ip)dip  =  0 
0  0 


a'  a' 

/  K  (ui,\p;y)  f  (ip)  dip  +  /  K  (oj,ip;y)g(ip)d\p  =  0 
0  Z1  0 

where 


Kn(a),\|>;y)  =  j  R0' 

K12(°J,^*y)  =  l'4 S * * 8m 

K21(w,^;y)  =  I'V 
K22(u^;y)  =  X’Hb* 

7 

It  can  be  shown  that 


+  l  Rjjj'  cos  m to  cos  imp 

cos  mio  sin  imp 

sin  mw  cos  imp 

sin  mw  sin  imp  . 


_m 

=  i^T  r  (2^£l)  r  (SiHii)  /r  (^)  /r  (2£^o 

Py 


(40) 

(41) 


(42) 


(43) 


Using  (43)  and  (34)  the  following  asymptotic  expansions  can  be  established. 

R  '  =  -(l-v-yv)_1(l/m)  +  0(l/m2)  ,  S  '  =  1/m2  +  0(l/m3) 

m  m 

U  '  =  r-(l-v)/(l+vy)] (l/m)  +  0(l/m2)  ,  T  '  =  1/m2  +  0(l/m3)  .  (44) 

m  m 

Substituting  from  (44)  into  (42)  and  summing  the  dominant  part  of  the  series 

/  4,8 

we  have. 


4Noble,  B.,  "The  Potential  and  Charge  Distribution  Near  the  Tip  of  a  Flat 

Angular  Sector,"  EM-135,  New  York  University,  NY. 

^Magnus,  W.  and  Oberhettinger,  F.,  SPECIAL  FUNCTIONS  OF  MATHEMATICAL  PHYSICS, 

Chelsea  Publishing  Co.,  1949. 

8Jolly,  L.B.W.,  SUMMATION  OF  SERIES,  Dover  Publications,  p.  126,  1961. 
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~  j  log  2  |  cos  to  -  cos  ijj|  +  regular  terms 


(45) 


K„„(o3,ijj)  =  —  log I  sin (~~^“) /sin (— ~^~)  I  +  regular  terms  (46) 

22  YJ  2  (1+vy)  1  2  2  1 

The  previous  expressions  show  the  kernels  have  logarithmic  singularities. 
Similar  analysis  can  be  carried  out  for  equations  (32)  and  (33)  for  mode  III 
(interchanging  R  1  and  U  '  S  '  and  T  ' ) . 

APPLICATION  OF  VARIATIONAL  PRINCIPLE 

Equations  (40)  and  (41)  are  identical  to  equations  (1)  and  (2)  and  all 

the  conditions  of  the  theorem  required  for  the  kernels  are  satisfied.  In 

this  section  we  shall  apply  the  theorem  to  obtain  approximate  eigenvalues  by 

assuming  approximate  trial  functions.  Without  causing  ambiguity  we  shall 

drop  the  asterisks  and  assume  the  following  trial  functions, 

1/2 

f(t)  =  (30  +  cos  t)cos(t/2) (cos  t  -  cos  a') 

1/2 

g(t)  =  [(1  -  cos  a')  +  2  cos  t] sin(t/2)/ (cos  t  -  cos  a')  (47) 

where 

n  2  r  ./-n1  +  U-V-vnORo*  +  log  sin2(a'/2) 

B  =  -  cos2 (a 72) - - -  (48) 

0  (l-v-vvi)Ro'  +  log  sin2(a'/2) 

The  above  trial  functions  are  the  first  approximations  to  the  integral  equa¬ 
tions  (40)  and  (41)  with  kernels  (42)  replaced  by  their  dominant  parts  given 
by  (44) .  The  method  of  obtaining  such  solutions  by  direct  computation  is 
illustrated  in  Reference  9. 


%oble,  B.,  Hussain,  M.A.,  and  Pu,  S0L0,  nApex  Singularities  for  Corner  Cracks 
Under  Opening,  Sliding,  and  Tearing  Modes, "  to  be  published  in  the  Proceed¬ 
ings  of  the  International  Conference  on  Fracture  Mechanics  in  Engineering 
Application,  Bangalore,  India,  1979. 
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Substituting  from  (47)  into  the  characteristic  equation  (3) ,  changing 

the  order  of  summation  and  integration  and  using  the  integral  representation 

_  7 

of  Legendre  functions,  we  obtain 


(49) 


where 


=  P  (cos  a') 
m  mv 


(50) 


For  different  values  of  v  and  a'  (the  complementary  angle  to  half  of 

the  vertex  angle  of  the  w'edge-shaped  crack)  we  get  approximate  values  of  y 

from  equation  (49)  for  both  modes  II  and  III. 

The  form  of  (49)  is  very  well  suited  for  Macsyma  evaluation  and  using 

eight  terms  for  summations  in  (50) ,  the  results  for  y  are  marked  by  x  in 

9 

Figure  3  where  the  solid  lines  are  results  obtained  by  another  method.  The 
results  by  both  methods  are  in  good  agreement. 

A  further  refinement  of  results  can  be  obtained  by  selecting 


f(t)  =  (A  +  B  cos  t) 


cos(t/2) 


(cos  t  -  cos  a*)l/2 


g(t)  =  (C  +  D  cos  t) 


sin(t/2) 


(51) 


(cos  t  -  cos  a')1^ 


7, 


Magnus,  W.  and  Oberhettinger,  F.,  SPECIAL  FUNCTIONS  OF  MATHEMATICAL  PHYSICS, 
Chelsea  Publishing  Co.,  1949. 
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and  formally  extending  the  variational  technique  to  a  characteristic  equation 
obtained  from  the  vanishing  of  the  determinant  of  a  four  by  four  system  (A, 

B,  C,  D  in  (51)  must  not  all  vanish).  The  results,  using  Macsyma  and  summing 
to  a  maximum  of  eight  terms  in  (50),  are  compared  in  the  following  table. 

They  also  are  shown  as  •  in  Figure  3. 


TABLE  II.  VALUES  OF  y  FOR  V  =  0.25 


Mode 

Half  Vertex 
Angle  a 

Variational  Method 

Using  (47)  Using  (51) 

Direct  Method 

Ref.  9 

O.lTT 

0,9333 

0.9616 

0.9582 

0.3tt 

0.6489 

0.6961 

0.6953 

II 

0.5tt 

0.4752 

0.5107 

0.5017 

0 . 7  TT 

0.3585 

0.3749 

0.3654 

0„9tt 

0.2213 

0.2469 

0.2137 

0.  ITT 

0.9739 

0o3tt 

0.7642 

0.8253 

0.8270 

III 

0.5tt 

0.4881 

0.5192 

0.5027 

0.6tt 

0.3829 

0.4034 

0.3914 

0. 7tt 

0.3021 

0.3015 

0o8tt 

0.2372 

0.2335 

9 

Noble,  B.,  Hussain,  M„A„,  and  Pu,  S.L.,  "Apex  Singularities  for  Corner  Cracks 
Under  Opening,  Sliding,  and  Tearing  Modes,"  to  be  published  in  the  Proceed¬ 
ings  of  the  International  Conference  on  Fracture  Mechanics  in  Engineering 
Application,  Bangalore,  India,  1979. 
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Even  the  rigorous  proof  of  the  extension  of  the  variational  technique 
from  the  two  by  two  system  using  trial  functions  (47)  to  the  four  by  four 
system  using  trial  functions  (51)  is  still  to  be  done,  the  results  thus 

obtained  are  in  good  agreement  with  results  achieved  by  using  trial  functions 

,  9 

(47)  or  other  methods „ 

CONCLUSIONS 

For  modes  II  and  III,  the  results  show  that  the  stress  singularities  are 
dominated  by  the  vertex  angle  as  well  as  the  elastic  constant  v  of  the  mate¬ 
rial.  The  results  further  indicate  that  when  the  apex  angle  is  greater  than 
180°,  the  stress  singularity  is  stronger  than  one  half  enhancing  the  tendency 
of  crack  front  to  straighten  out.  Similarly,  when  the  vertex  angle  is  less 
than  180°,  the  stress  singularity  is  less  severe  than  one  half  and,  again, 
this  will  tend  to  retard  the  growth  at  the  vertex  until  the  crack  front 
straightens  out. 

Macsyma  was  extensively  used  throughout  the  analysis,  especially  in  the 
generation  and  use  of  special  functions  such  as  Legendre  functions.  Gamma 
functions,  Bessel  functions  from  Share  directory,  in  the  summation  of  series, 
in  the  solution  of  linear  equations,  in  seeking  roots  of  determinants  of 
matrices,  in  the  plot  routine  and  in  the  creation  of  file  for  Batch  with 
Teco,  etc.  This  investigation  would  have  been  extremely  tedious  without 
Macsyma.  The  methods  as  well  as  results  in  the  full  entirety,  to  our  knowl¬ 
edge,  do  not  seem  to  have  appeared  in  literature. 

%oble,  B.,  Hussain,  M.A0,  and  Pu.  S.L.,  "Apex  Singularities  for  Corner  Cracks 
Under  Opening,  Sliding,  and  Tearing  Modes,"  to  be  published  in  the  Proceed¬ 
ings  of  International  Conf.  on  Fracture  Mechanics  in  Engr  Application, 
Bangalore,  India,  1979. 
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APPENDIX  1 


In  this  appendix  we  give  components  of  displacement  and  stress  in  terms 

of  Papkovitch  stress  functions  <J>,  \p,  w,  and  \0  A  subscript  to  a  stress 

function  means  the  partial  derivative  of  the  stress  function  with  respect 

to  the  variable  represented  by  the  subscript,  e<,g.,  4>r  =  9<})/3r,  $  = 

3(3<j)/3r)/3p.  The  variable  p  =  cos  9  is  used  in  the  place  of  0  in  both 

1/2 

Appendix  1  and  Appendix  2.  The  notation  p  =  sin  0  =  (1-p2)  is  also  used. 


The  First  Basic  Solution: 
2Gur  =  <J>r 
2Gu0  =  -  p<}>p/r 

2Gu<)>  =  V(rp} 


a  =  d> 

V  T 


cr  = 


o±  - 


rr 

r“2Cp2d)  -  pd>  +  rd)  ) 
PP  p  r 


0(f) 

rr(j) 


r  %4,/P  +  ^r  -  ^ 

-r~2  (^p<j)  +  P^/P2^ 

r'2P_1^r(}»  "  V 


t  0  =  r_  2P  C— _ +  4>J 


r0 


rp 


The  Second  Basic  Solution: 

2Gur  =  [r^»r  -  (3-4v)iJ;]p  cos  0 
2Gu0  =  -[p2^p  +  (3-4v)p^]cos  4> 

200^  =  ipQ  cos  <j>  +  (3-4v)iJ;  sin  <J> 

°r  =  tr^rr  "  2(1-v^r  +  2vr’lp,Jp]P  cos  ♦  + 


2v(rp^-1^ 


sin  (}) 
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Oq  =  r_1  [p2i|ipp  +  (l-2v)pi|>  +  (l-2v)npr]p  cos  <p  +  2v(rp)“1'JJ({)  sin  <j> 

Oq  =  Op)'1  +  (l-2v)rpz^r  -  (l-2v)pp2ijjp] cos  <p  +  2(l-v)  (rp)"1^  sin  <j> 

Td<p  =  OrP)_1  [p2^p(})  +  2(l-v)p^]cos  <p  -  (l-2v)r_1pjjp  sin  <p 
Tr<p  =  r_1tr^4,r  "  2(l-v)^]cos  4>  +  (l-2v)i|>r  sin  <p 
Tr0  =  ["P^rp  -  (l~2v)p^r  +  2(l-v)r-1p2^p] cos  <p 


The  Third  Basic  Solution: 


2Gur  =  [rw  -  (3-4v)io]p  sin  <J> 

2Gu0  =  -  [p2iop  +  (3-4v)pw] sin  <J) 

2Gu,  =  w.sin  d>  -  (3-4v)w  cos  d> 

<p  <p 

o r  =  [rco^  -  2(l-v)wi>  +  2vr_1pcop]p  sin  <J)  -  2v(rp)_1Lo^  cos  (f> 

Og  =  r_1  [P2(jJpp  +  (l-2v)pu  +  (l-2v)rwr]p  sin  <p  -  2v (rp) " 1  cos  <p 

=  r_1  [w^/p2  +  (l-2v)ro)r  -  (l-2v)pwp]p  sin  <J>  -  2  (1-v)  (rp)“  cos  <p 

Tqq  =  (-rpr1  [p2w  ^  +  2  (1-v)  pay]  sin  <p  +  (l-2v)r_1pwp  cos  <p 

T  .  =  r-1[rw  ,  -  2(l-v)co  ,]sin  <b  -  (l-2v)co  cos  <j> 
r<p  r<p  <J)  r 

Tr0  =  [-Plwrp  -  (l-2v)pwr  +  2  (1-v) r~ 1  p2top]  sin  <p 


The  Fourth  Basic  Solution: 

2Gur  =  [rXr  -  (3-4v)X]p 
2Gu0  =  [-pXp  +  (3-4v)X]p 
2Gu<p  =  pp-'X^ 

ar  =  rpXrr  "  2(1_v)PXr  "  2vr_1P2xp 

crQ  =  r_1pp2Xpp  +  (l-2v)pXr  +  r_1  (p2+2vp2-2)Xp 

a<p  =  r_1  [pP”2X(j)(j)  +  (l-2v)rpXr  -  (p2+2vp2)Xp] 
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Te<f>  -  r_1t-PVj>  +  C2-2v-p-)X4)] 

Tr(J>  =  -  2 (l-v)pA^] 

^r0  =  P["PXrp  +  ^ (1 -v) r” lpA  +  (l-2v)Xr] 
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APPENDIX  2 


In  the  following  solutions,  the  selection  of  cos  m<f>  or  sin  m<{>  depends 
on  the  geometry  of  the  problem.  The  sign  on  top  goes  with  the  trigonometric 
function  on  the  top  and  vice  versa. 


Solution  [A] : 

cos  m<}> 

2Gur  =  (y+l)ryp"V 

^  sin  m<}> 

2Guq  =  (ry/p)  [(y+l)pP™  -  (m+y+l)p"J] 


2Gu  =  ?  (ry/p)mP^+1 

cos  m<j) 

1  cos  m(f> 
a  =  yCy-H)^'^"1 
r  ^+^sin  m<J) 

a9  =  Oy  1/p2){m2-y-l-p2y(y+l)]p"|+1  +  p(m+y+l)P™} 


y 

sin  m<j) 


cos  m<j) 

I 

sin  nuj) 


a,  =  (ry  1/p2)  [(-m2+y+l)P™+1  -  p(m+y+l)P™] 


cos  m(j) 


cos  me}) 
sin  mej) 


t  =  +  (ry  1/p2)m[yPp"|+1  -  (m+y+l)P™] 


0(j) 

-  /  y-1 „m 

Tr(j)  =  +  O  /p)myP 

_y-i  <- 


m 

y-1 

sin  mef) 

y+1cos  mef) 
m 


sin  mej) 
cos  mef) 


sin  me}) 


TrQ  =  (r~  /p)y[(y+l)pP^+1  -  (m+y+l)p"J] 


Solution  [B] : 


cos  mej) 
sin  me}) 

cos  me}) 


2Gur  =  ry(y-3+4v)  [(-m+y+l)P^+1  -  (m+y+l)pP"J] 

sin  mef) 

2Guq  =  (ry/p) { (-m-4+4v) (-m+y+l)pPm+^  +  (m+y+1) [m+4-4v-p2 (y+4-4v) ]pm} 


cos  m<|) 
sin  me}) 
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2Gu,  *  *  (rP/p) (m+4-4v) [(-m+y+1) P™  .  -  (m+y+l)pPm] 

^  y+1 


sin  m<t> 


y  cos  mtj) 

•  m  m  cos  m<() 

a  -  r  [ (-m+y+1) (y-3y+2yv-2mv-2v)P  +  (m+y+1) (-yz+3y+2v)pP  ] 

r  y+1  y  sin  m<f> 


Og  ■  (ry  Vp2)  {  (-m+y+1)  [  (-yz-2yv- (3-2v)  (m+l)p2  +  (mz+(5-4v)m+4-4v)]P 


m 

y+i 


+  (m+y+l)p [-m2-5m+4mv-4+4v+p2 (yz+3y-3-2v) ] P  } 


0^  *  (ry  Vp2) {  (-m+y+1)  [-m2-5m+4mv-4+4v+p2  (l-2v)  (m+y+1) ]P™+^ 


cos  mtj> 

) 

^  sin  mtj> 


m 


+  (m+y+l)p [mz+5m-4mv+4-4v- (l-2v) (2y+l)p2] P  } 


cos  mtj) 


^  sin  mcj) 


T  =  ±  (r^  1/p2){(-m+y+l)p[m2+(5-4v)m+^(l-v)]pm 
0<P  y+1 


m  sin  mcf> 


+  (m+y+1) [p  (my+(3-2v)m+2(l-v)y+4(l-v))-m  -(5-4v)m-4(l-v)]P  } 

^  cos  m<i> 

_  y-1  mm  ^ltj) 

T  =  +  (r  /p) [m(y-2+2v)+2y(l-v)-2(l-v)] [(-m+y+l)P  -  (m+y+l)pP  ] 

r<^  ^  cos  mtj) 

xr6  =  (rP  1/p){(-m+y+l)p[-my+2(l-v) (m-y+l)]P^+1 


m 


+  (m+y+1) [my+2(l-v) (-m+y-1) +p2 (-yz+2-2v) ] P  } 


cos  mtj) 


y  sin  mtj) 


Solution  [C] : 

2Gur  =  ry(y-3+4v) (-P®+1  +  pP™) 


cos  mcf) 
sin  m (p 


u  m  cos  4 

2Guq  =  (rM/p){(-m+4-4v)pP^+1  +  [  (m-4+4v)  J-p2  (y+4-4v)  ]  P™} 


2Gu  =  +  (ry/p) (m-4+4v) (-Pm  +  pPm) 

y+i  y 


sin  ncf) 
cos  nc() 


a  ■  ry  1{ [-y2+y (3-2v)-2mv+2v]p"\  +  (y2-3y-2v)pPm) 


sin  mtj) 


cos  mtj) 


y+1 


y  sin  mtj) 
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P-1 


Of,  =  (r  /p2) { [p2 (yz+2vy- (3-2v) (m-1))  +  (-m2+(5-4v)m-4(l-v)) ]P 


m 


y+i 


m 


+  [p  (-yz-3y-3+2v)  +  m- (5-4v)m+4 (1-v) ]pP^} 


aA  =  1/p2)([p2 (l-2v) (m-y-1)  +  m2- (5-4v)m+4 (1-v) ]P™ 


cos  m<f> 
sin  me}) 


+  [pz (l-2v) (2y+l)  -  m2+(5-4v)m-4(l-v)]pPm} 


cos  me}) 


^  sin  me}) 


T0<J)  =  ±  (ry_1/p2)  [m2-(5-4v)m+4(l-v)]PP™+1 


+  [p2 (-my- (3-2v)m+2 (l-v)y+4 (1-v) ) -mz+ (5-4v)m-4 (1-v) ] P '} 


m. 


sin  me}) 


y 


cos  mej) 


Tr(J)  =  *  ^  1/P)  ^"my+ 


2  (l-v)m+2  (1-v)  (y-1)  ]  (-Py+1  +  pP™) 


sin  me}) 
cos  me}) 


Tr0  =  (ry"1/p){[-my+2(l-v)m+2(l-v)(y-l)]PP^+1 


+  [my-2(l-v)m-2(l-v) (y-l)-pz(-yz+2-2v)]P  } 


m  cos  me}) 


^  sin  me}) 


Solution  [D] : 

cos  me}) 

2Gu  =  (y-3+4v)rypPm 

r  r  u  •  , 

Hsin  me}) 


2Gu  =  (ry/p){(-m+y+l)pPm  +  [p2  (y+4-4v)  -  (y+1)  ]Pm} 
o  y+i 

sin  me}) 

2Gu  =  +  (mp/p)ryPm 

*  ycos  m<j) 

cos  me}) 

°r  =  ry  {2v(-m+y+l)P™  +  (y2-3y-2v)pP™> 

^  ^  sin  me}) 


cos  me}) 


^  sin  me}) 


aQ  =  (ry  X/p2) { (m-y-1) [l-(3-2v)p2]P^+1  +  [m2+y+l+p2 (-y2-3y-3+2v) ]pP™} 

y  y  sin  mi(i 

y  1  m  cos  m<£ 

a.  =  (r  '  /p2) {(-m+y+1) [1+p2 (-l+2v) ]P  +  [-m2-y-l+ (l-2v) (2y+l)p2]pPm} 

*  y+1  y  sin  m<t> 


cos  m<p 


21 


,  _  sin  m<J> 

T0  =  ±  (ry  /p2)m{(m-y-l)pP"‘  +  [y+2-p2(y+3-2v)]P™} 

^  M  ^  cos  m(jj 

,  sin  m<}) 

x  =  ±  (2-2v-y)mp(r  /p)P 

™  Mcos  m<j) 

Tr0  =  (rM"1/p)((-m+y+l}(y-2+2v)pPj;+1  +  [-y2+y  (l-2v)+2  (1-v) 

jn  COS  ITl{j) 

+  p2 (yz-2+2v) ] P  } 

y  sin  ratj) 
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Figure  1. 


A  flat  wedge-shaped  crack  under  three  different  inodes  and 
its  two-dimensional  counterparts 0 


■x 
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Figure  2 


A  spherical  coordinate  surrounding  the  apex  of  a  thin 
angular  sector  crack. 
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a /it 


Figure  3.  The  eigenvalue  y  as  a  function  of  cx/tt  for  v  =  0.25. 

(a)  Variational  results  indicated  by  x  and  •  for  mode  II. 
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Figure  3(b).  Variational  results  indicated  by  x  and  •  for  mode  III. 
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